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ABSTRACT 



We present detailed models of irradiated T Tauri disks including dust grain growth 
with power-law size distributions. The models assume complete mixing between dust 
and gas and solve for the vertical disk structure self-consistently including the heating 
effects of stellar irradiation as well as local viscous heating. For a given total dust mass, 
grain growth is found to decrease the vertical height of the surface where the optical 
depth to the stellar radiation becomes unit and thus the local irradiation heating, 
while increasing the disk emission at mm and sub-mm wavelengths. The resulting disk 
models are less geometrically thick than our previous models assuming interstellar 
medium dust, and agree better with observed spectral energy distributions and images 
of edge-on disks, like HK Tau/c and HH 30. The implications of models with grain 
growth for determining disk masses from long- wavelength emission are considered. 



Subject headings: Physical data and processes: accretion, accretion disks — stars: 
circumstellar matter, formation, pre-main sequence 



1. Introduction 



It has been clear for more than a decade that the absorbing surfaces or photospheres of T 
Tauri disks had to be "flared", i.e. curved upward away from the disk midplane, to explain the 
typical observed spectral energy distribution at infrared wavelengths (Kenyon & Hartmann 1987). 
Because of flaring, the irradiation of the disk by the central star tends to dominate the disk heating 
at large radii ( Calvet et al. 1992[ ). Since the amount of flaring depends upon the temperature 
structure, which in turn depends upon the flaring-dependent irradiation, the self-consistent 
structure in general must be solved iteratively (Kenyon & Hartmann 1987; D'Alessio 1996; 
D'Alessio et al. 1998, 1999a). The dust grains control the absorption and emission of radiation 
and therefore the knowledge of their properties help determine the disk structure. 
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In a previous paper (D'Alessio et al. 1999b = Paper II), we constructed models for T Tauri 
disks in which the vertical disk structure was self-consistently solved assuming well-mixed dust 
with properties characteristic of interstellar grains. We found that the resulting models were 
too geometrically thick, which is also the case of the two-layer flared disk models of Chiang 
& Goldreich (1997, 1999). Our models also exhibited too little mm- and sub-mm wavelength 
emission, in comparison with observations. These deficiencies can in principle be relieved or 
eliminated by models including grain growth and settling to the disk midplane, processes which are 



thought to be essential precursors to planet formation (e.g., Beckwith, Henning fc Nakagawa 2000 ). 
Even without settling, grain growth alone can help improve the comparison with observations, 
since for a given dust mass growth can increase the long-wavelength opacity at the same time it 
decreases the optical-near infrared opacity; the latter effect reduces the apparent disk thickness at 
short wavelengths, and can even reduce the true disk vertical scale height by reducing irradiation 
heating. 

It is generally thought that some grain growth must have occurred in T Tauri disks to explain 
the spectral indices of the long-wavelength disk emission in some systems (Beckwith & Sargent 
1991 = BS91; [Mannings fc Emerson 19941 ; |Pollack, et al. 1994| = P94). Unfortunately, the details 



of this process are poorly understood. We are therefore motivated to adopt a parameterized 
approach with few parameters to make an initial exploration of the effects of dust growth on T 
Tauri disk structure. 

In this paper we present sequences of detailed, physically self-consistent disk models with 
power-law distributions of dust particle sizes. We find that modest grain growth substantially 
improves comparison with spectral energy distributions (SEDs) and optical-near infrared imaging. 
Using both observed mm-wave fluxes and spectral slopes as constraints, we suggest limits on the 
maximum grain size in some systems. Our results emphasize the importance of optical depth 
effects in constraining disk properties, and suggest observational tests which might ultimately 
lead to more accurate disk masses. The presence of silicate emission may require differential 
distributions of small and large dust, i.e., settling of large grains to midplane with small grains 
(size < 10 /im) at larger heights. This issue will be addressed in a forthcoming paper of this series. 



2. Dust opacities 

Dust is the most important opacity source in T Tauri disk models, and controls the 
temperature distribution within the disk as well as the emissivity at long wavelengths. We adopt 
the model of dust for disks proposed by P94, with a few variations. We assume that the optical 
properties of olivine and orthopyroxene are represented by the "astronomical silicates" (Draine 



&: Lee 1984), i.e., crystalline olivine with dielectric functions given by Laor fc Draine (1993)| and 



Weingartner fc Draine (2000)[ However, following P94, the imaginary refractive index is taken to 



be constant for A > 800 fim, for consistency with laboratory measurements (Campbell & Ulrichs 



1969; Henning & Stognienko 1996). The optical constants of water ice are taken from Warren 
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1984. Refractive indices for troilite are taken from Begemann et al. (1994)1 (see penning et al 



1999) for 10 < A < 500 pm and from P94 for other wavelengths. Finally, the refractive indices of 
organics are taken from P94. The mass fractional abundances of these dust species respect to gas 
are Csii = 0.0034, £ tro = 0.000768, Cice = 0.0056, and ( org = 0.0041, and their bulk densities are 
p s ii = 3.3 , ptro = 4.83, pi C e = 0.92 and p org = 1.5 g cm -3 , respectively, consistent with P94. The 
sublimation temperature of each of these species is taken from P94. 



The efficiency factor Q a bs is calculated using a Mie scattering code ( Wiscombe 1979| , updated 



in 1996), assuming that the grains are spheres. Coagulation tends to produce grains which are 
not compact spheres, but rather aggregates with a porous structure and a bulk density that 
decreases with size ( Meakin fc Donn 1988| ). The shapes of these grains are affected by different 



kind of processes which may take place in the disk, such as turbulence, settling, radial drift, etc.; 



in turn, the optical properties of these aggregates depend on their complex shapes (e.g., Henning 



fc Stognienko 1996| ; [Beckwith, Henning fc Nakagawa 2000| ) . However, in the relatively dense 



environment of a circumstellar disk, there are frequent collisions between particles and a strong 
interaction with the surrounding gas, tending to produce erosion and compaction of grains (Miyake 
& Nakagawa 1993 = MN93). In any event, for simplicity and lacking any detailed description 
of the shapes and optical properties of grains in turbulent disks, we assume that the disk grains 
interact with light as compact segregated spheres. 

The size distribution of dust grains is taken to be a power law n(a) = uq a~ p , where a is the 
grain radius, uq is a normalization constant and the exponent p is a free parameter. Approximate 
justification for this is taken from the calculations of Weidenschilling (1997), who finds a wide 
range of dust sizes at any given time due to a combination of grain growth and shattering. We 
assume a minimum grain radius a m i n = 0.005 pm and different maximum size a max and exponent 
p. Table 1 summarizes representative values of the opacity and albedo at different wavelengths 
for various a max and p. In this paper we use the following notation: k u is the mass absorption 
coefficient (i.e., corresponding to "true absorptions"), Xv is the total extinction coefficient (i.e., 
absorptions plus scattering) and w u is the albedo (given by w u = 1 — k u /xu)- For comparison, 
we frequently refer to "ISM-dust", which corresponds to abundances and optical properties from 



Draine fc Lee (1984) , i.e., "astronomical silicates" and graphite, with a max =0.25 pm and p =3.5. 

Figure [l] shows the opacity as a function of wavelength for dust grain size distributions with 
p =2.5 and 3.5, and maximum size a max =10 pm, 1 mm and 10 cm, for temperatures T = 100 and 
300 K (i.e., with and without water ice, respectively). For a given value of p, when the maximum 
grain size increases, the opacity at small wavelengths (A <^ Q>max 

) decreases and the opacity at 
larger wavelengths increases. Since the dust to gas mass ratio is assumed to be constant and the 
mass in grains is dominated by the bigger grains if p < 4 (MN93), a larger a max implies a smaller 
mass in small grains which are the most efficient absorbing radiation in the visible and near-IR. 
At the same time, as a max increases, the transition between the geometrical and optical regimes 
occurs at longer wavelengths (~ 2-Ka max ), so that the opacity at A ~ a max increases. In Table 1 we 
show the Planck mean opacities characteristic of the stellar radiation k* p and x*pi both calculated 
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using the Planck function evaluated at a typical stellar temperature T* = 4000 K (see paper II). 
These values are similar to the monochromatic coefficient evaluated at A ~ 1 (im, and are used 
to calculate the disk heating by stellar radiation (see paper II for details). The table illustrates 
how these mean opacities also decrease with a max . The consequences of this behavior for the disk 
structure and emission properties will be discussed in §[|. 

Figure || shows the contribution of the different grains to the absorption coefficient for T = 
100 K, a max =1 /im and 1 mm and p = 2.5 and 3.5. In general, the relative contribution of each 
type of grain depends on the wavelength and size distribution. For a max = 1 mm, silicates are 
the most important contributor at A ~ 1 mm, and troilite becomes of comparable importance 
in the cm range. In the mid- and far-IR range, water ice dominates, and the second important 
contributor are organics. For & m ax — 1 A*^> silicates dominates the opacity in the mm and cm 
range, and troilite is not very important. In the mid- and far-IR water ice is dominant, with 
important contributions of silicates and organics. In particular, the organics tend to increase 
the opacity at A <J 6 fim, decreasing the contrast between the silicate band at 10 //m, and the 
adjacent continuum, respect to opacity calculations considering only silicates. The silicate bands 
at ~ 10 — 20 iim. almost disappear from the total absorption coefficient for a max = 1 mm (see also 
Figure [j]), even for temperatures higher than the ice sublimation temperature. We find a similar 
result for a max ^> 10 [xm ( in agreement with P94). 

Figure [|| (upper panel) shows the mass absorption coefficient at A = 1.3 mm as a function of 
maximum grain radius a max for p = 2.5 and 3.5, and for a temperature T < 100 K. The opacity is 
almost constant for a max < 10 /im, has a maximum around a max ~ 1 mm and decreases for larger 
Umax as k u oc a^P. for 3 < p < 4 and k u oc a^ ax for p < 3. This happens because the absorption 
coefficient at A ~ 1 mm is dominated by grains with a ~1 mm and the mass of those grains is 
maximum for a max ~ 1 mm and decreases for larger a max , with a slope that depends on p. Our 
results agree closely with those of MN93, even though we use a different dust composition. 

The straight line in Figure || (upper panel) is the opacity at 1.3 mm derived from the formula 
k v = 0.1(A/250 nm)~P , with (3 = 1 (BS91), commonly used to derive disk masses from mm-wave 
fluxes. At A =1.3 mm the value of this absorption coefficient is close to the maximum absorption 
coefficient, achieved when a max ~ 1 mm, but is generally much larger for other grain sizes. In 
addition, according to P94, composite, unfilled grains have in general smaller opacities at 1 mm 
than compact spheres of the same size. 

The frequency dependence of the dust opacity in the wavelength regime of interest, 
parameterized by /3 = dlog K u /d\og u, can in principle be used to constrain grain growth (BS91; 
MN93), depending on the chemical composition (P94). In Figure |3| (middle panel), we show the 
value of (3 in the mm-wave range, 

P _ log( K 0.769mm) ~ log(^1.3mm) 

log(1.3mm) — log(0. 769mm) ' 
as a function of maximum grain radius a max for the same cases plotted in the upper panel. For 
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small grains (a max < 10 /im), (3 = constant ~ 1.75, and approaches a maximum value for sizes 
a-max ~ 60 — 300 fim (corresponding to 2na/X ~ 1). The value of (3 decreases with increasing 
maximum size, and for a max > 10 cm, becomes independent of a max . Our results are qualitatively 
consistent with those of MN93. 

The detailed behavior of (3 for a max between 100 /im and 1 cm depends on the relative 
importance of the different kind of grains. Small changes in abundances and optical properties 
can change the value of (3 in these intermediate sizes. In particular, according to P94, composite 
grains with 50 % void volume have a less steep increase of (3 with a max , for grain sizes smaller 
than 1 mm. 

Figure || (lower panel) shows that the albedo oj v at A = 1.3 mm increases from to ~ 0.9 for 
a max between 10 /im and 1 cm (MN93). For larger values of a max , u> u ~ 0.9 for p=3.5 and oj u ~ 
0.5 for p =2.5. This large albedo means that scattered light contributions in the mm-wavelength 
range may not be neglected (MN93). 



3. Disk Models 

We solve self-consistently the complete set of vertical structure equations, including irradiation 
and viscous heating, resulting in detailed profiles of temperature and density with vertical height 
and disk radius. The assumptions, equations and method are described in Papers I and II. In 
particular, we assume that gas and dust are well mixed and in thermal balance everywhere. The 
input parameters to calculate the disk structure are the disk mass accretion rate M, the viscosity 
parameter a, and the central star mass M*, radius i?* and temperature T*. Whenever it is 
possible, we use observational constraints to reduce the parameter space that must be explored. 



To this end we have chosen typical parameters of T Tauri stars from Gullbring et al. (1998) : 
M = 3 x lO^Moyr- 1 , M* = 0.5 M , R* = 2 R and T* = 4000 K for the fiducial model. We 
assume the disk is in steady state accretion, i.e. the accretion rate through the outer disk is the 
same as in the inner disk, at a rate inferred from measurement of the ultraviolet and optical 
continuum excess emission produced as material accretes via a shock at the stellar surface ( Calvet| 



and Gullbring 1998|) . The viscosity parameter a is changed appropriately to fix the disk mass. 



The adopted values of a are between 0.005 and 0.036, corresponding to disk masses between 0.1 
and 0.02 M (for M = 3 x 10~ 8 A^yr" 1 and R d = 100 AU). The disk outer radius is taken to be 
Rd = 50, 100 and 200 AU, and the inner radius is Rhole = 3 i?*. 

Given the disk structure, the emergent intensity is calculated by integrating the transfer 
equation, 

I v = I b J> e" T - + f " S u (t Vii ) dt Vti , (2) 



o 



along rays that pierce the disk. Here, S u is the source function, t v ^ is the total optical depth along 
a particular ray (with an inclination angle i respect to the symmetry axis of the disk), t„ $ is the 
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variable optical depth along the ray, and I„ 9 is the background radiation, i.e., I^ 9 = B u (Tb g ) with 
Tj, g = 2.73 K. The optical depth is calculated including true absorptions and scattering events. 

Since the dust albedo of large grains is so large in the mm range (Figure ^ and Table 1, it is 
important to include the contribution of the emissivity in scattered light. Assuming coherent and 
isotropic scattering, the source function is taken to be QMihalas 1978| ) 

S v = (1 — w v )B v + u u J v , (3) 

where B u is the Planck function and J v is the mean intensity of the radiation field. We approximate 
the J v by the mean intensity of a vertically isothermal slab (see MN93), with the Planck function 
evaluated at the local temperature, 

T x on Ji i exp t~ v /3 ( 1 ~ Ml + ex P [v /3 ( 1 ~ u v) (%Q ~ T ".o)] 1 f A \ 

Mtufi) = B v {t v>0 )<l + ■ — > , (4) 

I exp [— v / 3(l - u>„) r^oJivl - ^ — 1) — (vl - w„ + 1) J 

where r^o is the total optical depth in the vertical direction (i = 0°) and t^o is the variable optical 
depth in the same direction. 

In an attempt to simulate the observational procedure of taking on-off measurements, the 
background intensity I^ g due to the cosmic microwave background is subtracted from the disk 
emergent intensity given by equation (§). The observed flux is calculated by integrating the 
intensity, emerging in a direction defined by the line of sight, over the solid angle subtended by 
the disk as seen by the observer. For optically thin regions, the emergent flux calculated including 
scattering is similar to the flux calculated assuming the dust is purely absorbing. For optically 
thick regions, with r^o >> 1/(1 — u> u ), the flux for a pole-on disk is smaller by a factor roughly 
1 — u) v j\\J\ — u> u + 1] [ a/3(1 — u u ) + 1] than the purely absorbing case, consistent with the discussion 
in MN93. For intermediate optical depths (e.g. t u> q ~ 1), the behavior of the emergent flux cannot 
be described by analytic expressions. 



4. Results 

4.1. Disk "photosphere" and hidden central stars 

Figure [l] shows that grain growth affects the opacity at both small and large wavelengths. For 
a fixed mass of dust, increasing a max tends to increase the mm- wave opacity (for a max mm) at 
the same time as it decreases the opacity at the characteristic wavelength of the stellar radiation, 
A ~ 1 /iin. This has important effects on the disk temperature distribution and thus the disk 
structure. 

Figures ||] and || show the midplane temperature T c , surface density S and the height of 
the "irradiation surface" function of radius, and the SEDs of pole-on (i = 0°) disks, 

for different grain size distributions and the same disk mass ~ 0.046 M Q , accretion rate 
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M = 3 x 10~ 8 M© yr -1 and the typical central star properties mentioned in §|3[ The height z s is 
defined to be where the mean optical depth to the stellar radiation is unity (using the mean opacity 
Xp, listed in Table 1) and represents where most of the stellar radiative energy is deposited (see 
paper II). The height of this surface decreases when the grain sizes increases. Thus, the fraction 
of stellar radiation intercepted by the disk decreases, resulting in a colder disk than what would 
be expected for smaller grains. The mass surface density of the outer disk is S ~ 1/R for R > 
10 AU, thus the disk mass is dominated by the contribution of the outer radii. Since the present 
models have all the same disk mass, they also have the same surface density for R ;> 10 AU. 

In Paper II we showed that if the disk has ISM-dust well-mixed with the gas, then a large 
fraction of the CTTS should have their central star extinguished with Ay > 30 mag (30 % 
of all objects with an outer disk radius Rd = 100 AU). Also, the two-layer models of Chiang 
& Goldreich (1997,1998) implies a similar large fraction of highly-extinguished central stars. 
As discussed in Paper II, this fraction is too high in comparison with estimates from current 
observational surveys in the Taurus molecular cloud complex. The decrease in the optical and 
near-infrared dust opacity with increasing a max results in (optically and geometrically) thinner 
disks, and thus the fraction of highly-extinguished stars predicted for a random distribution 
of viewing inclinations decreases. Figure ^ shows the cosine of the critical angle i c , such that 
cos i c = /x c = za v {R)/[za v {R) 2 + .R 2 ] 1 / 2 , where za v {R) is the height at which the the extinction 
coefficient is Ay = 30. Assuming a random distribution of disk rotation axes to the line of sight, 
for a max = 1 mm and p = 3.5, 20 % of the T Tauri stars with disks with a typical radius Rd ~ 100 
AU, should be extincted by their disks with an Ay > 30 mag, while for a max = 10 cm and p = 2.5 
only ~ 6 % of T Tauri should have Ay > 30 mag. Because the fraction of such disk-obscured 
sources in Taurus is estimated to be not more than about 15 % (Paper II), it is clear that the 
grain-growth models are in much better agreement with the current observational constraints. 

The SEDs of the models show that the larger the fraction of big grains, the smaller the mid- 
and far-IR fluxes. On the other hand, the mm fluxes increases with a max until a max ~ 1 mm, and 
decreases with a max for larger grains, reflecting the behavior of the opacity shown in Figure ^. In 
the next section we show how the SED of a disk model with a larger fraction of big grains than 
the ISM dust compare to observations. 

4.2. Median SED 

Figure shows the spectral energy distributions (SEDs) for two disk models, both with 
M d = 0.046 M© (M = 3 x 10~ 8 M Q yr-\ a = 0.01, R d = 100 AU) and a central star with 
M* = 0.5 M , i?* = 2 R and T* =4000 K (L# = 0.9 L Q ), and an inclination angle cos i = 0.65 0. 



4 We use cos i — 0.65 instead of 0.5 for the average inclination to account for the occultation in edge-on systems, 
which would drop out of the observational samples due to weak optical and infrared emission. 
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We adopt the distance to Taurus of 140 pc ( [Kenyon, Dobrzycka and Hartmann 1994 ). One of the 
models has ISM-dust and the second model has dust with the abundances and optical properties 
described in §|2], and a max = 1 mm, p =3.5. The main difference between both models are that 
the ISM-dust disk model has larger far-IR and smaller mm-wave fluxes than the grain growth 
model. Both models are compared to the median SED for the T Tauri stars in the Taurus-Auriga 
molecular cloud ( from Kenyon fc Hartmann 1995 , c.f., Paper II). The SED for the a max = 1 mm 
disk model is much closer to the observed median spectrum. A range of values of a max and p = 
2.5 - 3.5 produce SEDs consistent with the quartiles (c.f. Figures § and ||). Also, disks must have 
different radii, inclination angles, environment, etc., which would contribute to the scatter in the 
observed SEDs (c.f., paper II). 



4.3. Scattered light images 

The decrease in opacity at small wavelengths when the fraction of large grains increases 
also affects images in stellar scattered light. Figure |8] shows images in scattered light of edge-on 
disk models (i = 90°), with the typical properties described in §|| and a viscosity parameter 
a = 0.01. These disk models have Rd = 100 AU and Md = 0.046 M Q , and again are assumed to 
lie at a distance of 140 pc. The images are calculated at A = 0.814 fj,m, in the single scattering 
approximation, assuming a Henyey-Greenstein phase function, with an asymmetry parameter 
g = 0.65. The emergent intensity distribution has been convolved with the HST PSF (for the 
filter F814W/PC1). The aspect ratio of the image, i.e. the width of the central dark lane or 
equivalently the distance between emission maxima, is very sensitive to a max for a given p. In this 
case, the larger the value of a max the smaller the width of the dark lane, which almost disappears 
for p = 2.5, a max =10 cm. This occurs because the opacity at A ~ 0.8 pm decreases when a max 
increases and p decreases (c.f and Table 1), and both the height of the surface where the 
stellar radiation is attenuated (c.f §^4]) and the extinction of the scattered light at a given height 
produced by the outer disk decrease. For the same reason, the maximum intensity of the images 
increases with a max . The width of the dark lane also depends on the outer disk radius. The height 
where most of the stellar radiation is deposited (z s , c.f. Figures Q and |5|) evaluated at the the disk 
radius can be used to make a rough prediction of the distance between the emission peaks. It is 
not very sensitive to the disk inclination, as long as the the disk occults the central star. 

In Paper II we showed that the near-IR scattered light image predicted for the well-mixed 
ISM-dust models produces an aspect ratio ~ 2 times larger than is consistent with Hubble 
Space Telescope images of HK Tau/c and HH 30. The results shown in figure || show that grain 
growth will improve the comparison between theory and observation. Here we consider a detailed 
comparison between the the image of HK Tau/c and a disk model. 

Figure || shows the image of a disk model with M = 3 x 10 -8 , a = 0.01, T* = 3500 K, 
R* = 1.3 R (£* = 0.22 L Q ), M* = 0.5 M , R d = 110 AU, inclination angle i = 89.8° and position 
angle P.A.= 40° . The disk model dust is characterized by a max = 1 m and p = 3.5, and the total 
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opacity and albedo are given in Table 1. 



The model image is compared to the image of HK Tau/c ( Stapclfeldt et al. 1998| ), with the 



same contour levels in both images as in Figure ||[ The observed image and the model show 
reasonable agreement, as for instance the distance between maxima, the aspect ratio and the 
maximum brightness. The flux at A =1.3 mm predicted by this model is F\$ mm = 8 mJy, 
consistently smaller than the flux for the unresolved HK Tau binary system reported by [Beckwith 



et al. (1990)| (= BSCG), F L3 mm = 41 mJy. The mass of this disk model is M d = 0.065 M . We 
find that growth to a very large value of a max yields a dramatically reduced mm-wave opacity, and 
thus a very low mm-wave flux for this large disk mass (see §5), while producing an optical image 
in good agreement with observations. 

On the other hand, dust grains with a max = 1 mm and p = 2.5 have, approximately, the same 
opacity and albedo at A = 0.8 fim, but a much larger opacity at 1.3 mm (c.f., Table 1). We made 
a disk model with the same mass and central star properties than the one shown in Figure || and 
we found it has a similar scattered light image, but a much larger flux at 1.3 mm, Fi^ mm ~ 30 
mJy. However, this flux seems too high, considering that the total flux should have contributions 
from both disks in the binary, and that the less inclined disk of the primary is expected to have a 
larger contribution. 

The model shown in Figure |9] is not unique, being under-constrained by observations. 
Nonetheless, this example illustrates how disk masses inferred from optical imaging (or even 
mm-wave fluxes) may be substantially under-estimated. Without additional observational 
constraints on dust properties, detailed conclusions are difficult to derive. 

Moreover, the mass of the present model is several hundred times larger than the mass inferred 
by Stapelfeldt et al. (1998) (i.e., their models C and D have masses M d ~ 5.9x 10~ 5 -1.2x 10~ 4 M Q ). 
The discrepancy is due to the difference in dust opacity. Stapelfeldt et al. used an opacity typical 
of ISM dust grains, which is much larger at A = 0.8 fim than the opacity corresponding to a larger 
value of a max (see Table 1). This implies that they require a much smaller disk mass than we do 
to reproduce the width of the dark lane of the image and its aspect ratio. Assuming that the disk 
is optically thin and the flux at 1.3 mm scales with disk mass and the absorption coefficient, using 
the values given in Table 1 we infer F\^mm ~ 0.01 mJy for their model D. Thus, mm observations 
resolving the binary system would be very important to define whether this is a very low-mass 
disk with ISM type of dust or a moderate mass disk with a much smaller fraction of small grains 
than the ISM dust, as we suggest here. 



5. Disk mm-wave fluxes and masses 

Much of what is known or estimated about T Tauri disks comes from analyses of mm- and 
sub-mm wavelength emission. The first observational arguments in support of grain growth were 
made to explain observed spectral indices in this wavelength regime (BSCG, BS91). Our previous 
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models (Paper II) with ISM-dust could not explain the observed mm-wave fluxes. In this section 
we discuss our detailed model results for disks with larger grains in this wavelength region. 

In principle, disk masses and aspects of grain properties can be derived from long- wavelength 
fluxes and spectral slopes. The analysis is particularly easy when the disk is optically thin. In this 
case the observed flux is proportional to the disk mass, absorption coefficient at the wavelength 
of observation, and the source function at this wavelength. At long wavelengths the source 
function is expected to have a Rayleigh-Jeans wavelength or frequency dependence, so that the 
flux is then F u oc k (T}/c 2 v 2 k u , where k is the Boltzmann constant, c is the speed of light, 
and (T) is a mass- weighted disk mean temperature. The spectral index dlog F u /d\og v is the 
product of the intrinsic dependence of the dust opacity times the Rayleigh-Jeans distribution. 
Unfortunately, optical depth effects cannot be ignored in general and the source function may 
not have a Rayleigh-Jeans form if the dust temperature is too low (BS91), so in general detailed 
models must be used to compare with observations. 

To make this comparison we constructed model sequences of fixed disk mass with differing 
dust properties. The stellar luminosity, stellar mass, and radius are fixed as before at 0.9 L©, 0.5 
M© and 2 R (c.f., §||). We assumed a fixed mass accretion rate M = 3x 10~ 8 Mq yr _1 and a fixed 
outer radius of Rd = 100 AU in our basic sequence. Changing the disk mass in our models then 
requires us to change the a parameter accordingly. Since the viscous heating in the outer disk in 
the basic sequence is small, so that the temperature distribution is controlled by irradiation, this 
is roughly equivalent to adopting a constant viscosity parameter and varying M instead to change 
the disk mass. The disk masses in each sequence are = 0.023 Mq, 0.046 Mq, and 0.092 Mq. 
We assume initially that the disks are viewed pole-on, and discuss the effects of changing the 
inclination angle later. 

Figure |l0] shows the 1.3 mm flux vs. n for several sequences of disk models of fixed mass 
(connected points) The spectral index n is given by 

n = log(^0.769mm) ~ logt^Umm) , g s 

log(1.3mm) — log(0. 769mm) 

Each point along a given sequence corresponds to a different value of a max (between 10 fim and 1 
cm); the two panels correspond to p = 3.5 (upper) and p = 2.5 (lower). 

We discuss the model sequences beginning on the right-hand side of the figure, where the 
models have small a max , and moving counter-clockwise along the curves as a max increases. For 
small a max , the spectral index is relatively large, reflecting the intrinsic wavelength dependence of 
the dust opacity, though the actual value of n departs from the Rayleigh-Jeans limit 2 + (3 ~ 3.75 
because this limit is not applicable in this wavelength range for the low temperatures (~ 10 — 20 K) 
characteristic of the outer disk. The flux levels scale nearly linearly with the mass, as expected in 
the optically-thin case; only the inner, high-surface density regions are optically thick, and these 
regions do not contain much mass. Proceeding to larger drnaxi 

the spectral index decreases, both 
because larger regions of the disk become optically thick due to the increase in opacity (Figure |3|) 
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and because of the decrease in j3, which affects the frequency dependence of the emission from 
the optically thin regions (cf. Figure |3|). As a consequence of the grain albedo, the total optical 
depth of the disk has to be tq ;> 1/(1 — w) 1 / 2 for the emergent flux to approach the optically thick 
limit for pure thermal emission. For the models plotted in Figure |l0|, this limit is reached only 
by the most massive disk (Md ~ 0.1 M ) when a max = 1 mm. Finally, for much larger a max , 
the disk models become mostly optically thin, and so the fluxes and spectral indexes much more 
closely reflect the dust opacity and disk mass. The p = 2.5 models show much more variation in n 
at large dust sizes because (3 changes more rapidly with increasing a max than in the p = 3.5 case 
(Figure [|). For reference we also show the flux and spectral index of disk models with the same 
masses calculated assuming they contain ISM-dust (large open triangles). In this case the outer 
disks are optically thin, and the opacity is too small to account for the observed fluxes. 

To compare with our models we show observations of Taurus stars taken from BSCG and 



BS91 as stars in Figure |10|. The observed spectral index is calculated as n = j3 p + 2, where 



taken from Beckwith & Sargent (1991), is the result of power law fits to the spectral energy 



distributions between 1.3 mm and either 0.624 mm or 0.769 mm. We ignore the effects of the 
different wavelength limits on the spectral index; this has little effect on the comparison, given the 
typical errors shown by the error bar in Figure (BS91). Fi gure id shows that the disk models 



generally span the range of observed fluxes and spectral indexes larger or similar to n ~ 2. 



Figure |10| also shows the flux and spectral index of flat irradiated disk models, with a max = 1 
mm. Here we have assumed that the disks are geometrically flat, which corresponds to complete 
dust settling to the disk midplane. These disks are much colder at a given radius than those in 
which dust and gas are well mixed up to several scale heights (flared disks), since they intercept a 
smaller fraction of the stellar radiation (e.g., Hartmann 1998| ). The emission of these flat disks is 



dominated by optically thick and cold regions and both the flux and the spectral index are smaller 
(n < 2) than in the flared models. Thus, the CTTS with n < 2 and a larger mm flux than the flat 
disk models could be intermediate cases, in which just a fraction of dust has been settled towards 
the midplane. 

Panel (a) in figure |ll| shows the effect of changing the outer disk radius. We have selected the 
parameters (i.e., M and a) of the model with intermediate mass in Figure |l0[ The models have 
Rd = 50 AU (dashed line), Rd = 100 AU (solid line) and Rd = 200 AU (dotted line), corresponding 
to masses Md ~ 0.02, 0.05 and 0.092 M©, respectively (similar to the disk masses represented 

I — K 1/2 

in Figure |10|). In general, the flux scales as ~ R d , as expected for an optically thin disk with 
I v ~ B v t u oc i?~ 3 / 2 . However, for a max ~ 1 mm and Rd = 50 AU, a larger fraction of the disk 
becomes optically thick which is reflected by n — > 2. 



Panel (b) in figure 11 shows the effect of changing inclination angle. For the same model with 
intermediate mass and a radius Rd = 100 AU (M^ ~ 0.046 M ), we change the inclination angle 
cos i = 0.25 (dashed), 0.5 (solid) and 0.9 (dotted). For a max < 150 pm, the fluxes are similar 
reflecting that the disk is mainly optically thin at A = 1.3 mm, but the spectral index tends to 
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decrease with the inclination of the disk, because the relative contribution of the optically thick 
inner regions changes at both wavelengths. For a max > 1 mm, the flux roughly scales with cos i, 
and the spectral index is almost constant n ~ 2.4. 



Panel (c) in figure 11 shows the flux vs spectral index of pole-on disk models with different 
mass accretion rates and masses, assuming no irradiation by the central star, i.e. the disk heating 
is entirely due to viscous dissipation (Paper II). Each connected sequence corresponds to disk 
masses of = 0.023, 0.046 and 0.092 M and each point along a sequence corresponds to a 
different mass accretion rate (roughly differing by a factor 3). We take a max = 1 mm and p = 3.5 
for all models. The flux increases with temperature (or M) and the spectral index is sensitive to 
both the disk temperature and mass. In these models the outer disk is very cold (T ~ 5 — 15 
K)[] and the Rayleigh-Jeans limit is not valid. Unfortunately, since the optical depth is close to 
unity, one cannot use neither the optically thin nor the optically thick limits of equation (0) to 
approximate the mm-wave fluxes emerging from the disk by simple expressions. In addition, the 
disk temperatures are comparable to the cosmic background temperature and the optical depth of 
the disk is large, so the correction introduced when we substract B u (Tb g ) to the emergent intensity 
of the disk (c.f. §|3|) becomes very important. Besides this complications in describing analytically 



the results shown in Figure [Ll], the conclusions are simple. Larger mass accretion rates for these 
pure viscous disk models result in larger mm-wave fluxes. On the other hand, the smaller the disk 
mass, the smaller the flux and the larger the spectral index, because the smaller the optical depth. 
These models have 1 n 2.6, but require large mass accretion rates and small inclination angles 
to have mm-wave fluxes consistent with the observations. 

Panel (d) shows the same models as in panel (c), but including irradiation from the central 
star. All the models with different mass accretion rates at a given mass have almost the same flux 
at 1.3 mm and the same spectral index, because the stellar irradiation is the dominant heating 
source of the outer disk (R ^ 20 AU) . 



6. Discussion 

We have shown that self-consistent disk models with power-law distributions of grains, with 
p ~ 2.5 — 3.5, and maximum sizes a max ~ 1 mm, can reproduce many general features of observed 
T Tauri disks, including SEDs, optical and near-infrared images, limits on the frequency of 
disk-obscured systems, and mm-wavelength fluxes and spectral indexes. In contrast, our results of 
Paper II showed that ISM dust properties were not in good agreement with observations. Our 
findings reinforce previous suggestions (e.g., BS91, MN93) for grain growth in T Tauri disks with 
typical ages of 1 Myr. We make these suggestions stronger, however, since we aim to fulfill more 



5 We include the disk heating by cosmic rays and radioactive decay of 26 Al, but the disk is so dense that it cools 
efficiently, and can reach temperatures lower than 10 K. 
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observational constraints. 

The inference of substantial grain growth implies that there must be significant uncertainties 
in deriving disk properties such as masses. If dust agglomeration has proceeded from typical 
interstellar sizes of 0.1 /iin to a 1 mm or larger, there is no particular reason to assume that 
growth has proceeded to precisely the same value of a max in all systems. While it may be possible 
to argue that a max is at least ~ 1 mm on the basis of observed spectral indexes (Figure 8), it is 
difficult to eliminate the possibility of growth beyond this size from n alone because (3 approaches 
a constant for larger a max (Figure 3). This uncertainty in dust size distributions (even assuming 
p can be constrained to lie between ~ 2.5 — 3.5) results in a corresponding uncertainty in disk 
masses and other properties. 

Limiting dust opacity variations by constraining a max to lie between 100 //m and 1 cm, our 
disk models suggest that disk masses between ~ 10~ 2 Mq and ~ 10 are consistent with the 

data. The upper end of this mass range encompasses disks subject to gravitational instability. 
Specifically, the Toomre parameter Q (e.g., D'Alessio et al. 1999b) at Rd = 100 AU is Q ~ 1 
for the Md = 0.05M Q models and Q ~ 0.5 in the Md = O.1M models, suggesting that the 
former model is close to instability and the latter is strongly unstable. Thus, the possibility of 
gravitational instabilities in at least some T Tauri disks cannot be dismissed. In this regard, 
models shown in panels (c) and (d) of Figure |ll| indicate that observed high mm fluxes and 
low spectral indices of some of the stars could be explained by high mass disks with high mass 
accretion rates, as long as their outer disks were not heated by irradiation. Interestingly, some 
of the stars with n < 2 (CW Tau, DR Tau and DL Tau) are high mass accretion rate stars with 
M ~ 2 x 10~ 7 — 8 x 10 -6 M yr _1 ([Hartigan, Edwards fc Ghandour 1995[ ). We may speculate 



that one possible interpretation is that these disks may be gravitationally unstable and the outer 
disks are shadowed from the star by density perturbations produced by the development of the 
instability, as models suggest (e.g., Pickett et al. 2000, Nelson et al. 2000). If one assumes 
that the disk masses are always smaller than the gravitational stability limit, then the high mm 
fluxes shown in Figure cannot be understood if grains have grown beyond about 1 cm in size. 
However, growth to larger sizes could be accommodated in the inner disk; these are optically thick 
regions and the resulting emission is insensitive to their total mass. 

The models predict that less-luminous disks are more likely to be optically thin and thus 
exhibit spectral indexes which reflect the intrinsic wavelength dependence of the dust opacity 
much more closely. There are many nearby T Tauri stars whose 1.3 mm fluxes lie in the range 
of ~ 15-30 mJy (psterloh fc Beckwith 19951) for which mm-wave spectral indexes have not 



been published. Observations of spectral indexes for these fainter objects could more clearly 
demonstrate particle growth. However, because the spectral index tends to a limiting value when 
a max becomes large, it will be difficult to tell whether low mm-wave emission is due to low intrinsic 
disk mass or particle growth, unless mid- and far-IR fluxes and scattered light images are used to 
disentangle disk mass and grain sizes. 
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In principle, geometrically-thinner disks could be the result of dust settling alone. However, 
our analysis of the mm-wave fluxes and spectral indexes indicates that grain growth is essential. 
Because the growth of particles reduces the opacity at short wavelengths, it is difficult to 
distinguish between the effects of settling and grain growth on observed SEDs and images. The 
case of HK Tau/c illustrates how one may arrive at very different conclusions about disk masses 
depending upon the dust properties assumed. 

Of course, grain growth should be accompanied by settling; larger particles tend to fall more 
rapidly toward the disk midplane (e.g., Weidenschilling Sz Cuzzi 1993 , Weidenschilling 1997| ). 



Miyake fc Nakagwa (1995)| argue that settling has to occur in some T Tauri disks on the basis 
of their infrared SEDs. Two of our results also support the idea of settling. First, our models 
generally cannot explain the T Tauri stars with moderately low mm-wave fluxes and small values of 



n (e.g., Figure 10). Flat disks produce lower mm-wave fluxes and smaller spectral indexes, because 
the reduced irradiation fluxes reduce the outer disk temperatures. Our detailed calculations show 
that geometrically- flat disks result in too small values of flux and n, but the results shown in 



Figure Figure 10 suggest that disks with intermediate thicknesses (i.e., some settling of dust) could 



explain the T Tauri disks found in the lower-left hand region of the flux-spectral index diagram. 

Second, our models with a max ~ 1 mm and p = 2.5 — 3.5 have no silicate feature, in contrast 
with observations of many T Tauri stars ( |Cohen and Witteborn (1985) , Natta, Meyer &: Beckwith 



2000[) . In the context of our disk models, in which the relative grain size distribution is independent 
of height and radius, putting more of the total dust mass into large grains (grain growth) reduces 
the geometrical disk thickness at large radii, in better agreement with observations, but results 
in too few small grains at 1 AU (T ~ 100-300 K) to produce 10 fi silicate emission. We could 
improve the agreement with observations by having more small grains at 1 AU than at large radii. 
If the grain size distribution is independent of height, this would imply that there has been less 
grain growth at small radii than at large disk radii. It seems much more likely that there has been 
faster grain growth at small radii, but that the large grains have also settled more toward the 
midplane in the inner disk. This settling could result in a much larger fraction of small grains in 
the upper atmospheric layers at 1 AU, which is the region where the silicate emission is formed 
dCalvet et al. 1992| ; |Chiang & Goldreich 1997|) . 



It is therefore clear that improved disk models for T Tauri stars will have to incorporate both 
grain growth and settling in a physically self-consistent manner, accounting for the radial and 
vertical dependences of both process. We will develop such models in a forthcoming paper. 
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Fig. 1. — Mass absorption coefficient k v vs wavelength for a grain size distribution characterized 
by P = 3.5 (upper panels) and 2.5 (lower panels). Each curve corresponds to a different grain 
maximum size: a max = 1 /ira (dashed line), 1 mm (solid line), 10 cm (dotted line). The absorption 
coefficient is shown for T = 100 K (i.e, all the different kind of grains are present) and T = 300 K 
(i.e., water ice has been sublimated). Each panel is labeled by p and T. 
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Fig. 2. — Mass absorption coefficient k u vs wavelength for a grain size distribution characterized 
by P = 3.5 (upper panels) and p = 2.5 (lower panels), with a max = 1 //m (left panels) and 1 mm 
(right panels), for a temperature T = 100 K. Each curve corresponds to a different kind of grain: 
silicates (light solid line), water ice (dotted line), organics (dashed line), troilite (dot-dashed line). 
The total mass absorption coefficient is represented with a solid thick line. Each panel is labeled 
by p and a rnax . 
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Fig. 3. — Upper Panel: Mass absorption coefficient at A = 1.3 mm as a function of the maximum 
grain radius a max . The solid line corresponds to p = 3.5 and the dotted line to p = 2.5. The 
horizontal dashed line represents the frequently adopted opacity at 1.3 mm from BS91. Middle 
panel: (3 = dlog(K I/ )/dlog(i^), calculated between A = 0.769 and 1.3 mm for the same cases as the 
upper panel. The horizontal dashed line is (3 =1. Lower panel: albedo at 1.3 mm for the same 
cases shown in the upper panel. 




Fig. 4. — Disk structure and SED for different maximum grain size and exponent p=3.5. Each 
curve corresponds to a different value of maximum grain radius, a max = 10 fim (dashed line), 1 
mm (solid line), 1 cm (dot-dashed line) and 10 cm (dotted line). Each panel represents the radial 
distribution of: (a) Midplane temperature; (b) Mass surface density; (c) Height of the irradiation 
surface divided by radius. In panel (d) the resulting SEDs are shown, assuming the disks are 
oriented pole-on. The SED of the central star is also plotted (long-short dashed line). 
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Fig. 5. — Disk structure and SED for different maximum grain size and exponent p=2.5. See 
caption of Figure |j. 
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Fig. 6. — Cosine of critical inclination angle for disks models with = 0.046 M© (M = 
3 x 10~ 8 M Q yr _1 , a = 0.01). The surface is defined so the extinction of the central star produced 
by the disk is Ay = 30, for dust abundances and optical properties described in §|2| and different 
distributions of grain sizes: a max = 1 mm (solid line), 1 cm (dotted line) and 10 cm (dashed line). 
The upper panel corresponds to p = 3.5 and the lower panel, to p = 2.5. For comparison we also 
show [i c for ISM-dust, with ingredients and optical properties from Draine Sz Lee (1984)| , p = 3.5, 
o-max = 0.25 [ira (dot-dashed line). 
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Fig. 7.— SED of the model with M d = 0.046 M , M = 3 x 10~ 8 M yr _1 , cosi = 0.65, R d = 100 
AU, for ISM-dust (dotted line) and for the dust model adopted in this paper, with p = 3.5 and 
Umax = 1 mm (solid line). The points represent the median observed SED and the errorbars are 
the quartiles. The observed fluxes are normalized at A = 1.6 fim. For the mm wavelengths we also 
show the median fluxes without normalization (crosses) (normalization is a reasonable procedure 
for optically thick emission, but since the mm flux could have an important contribution from 
optically thin region, the true mm median would be something between the normalized and the 
non normalized fluxes). The contribution of the UV and optical excess emission produced in an 
accretion shock is not included in the model SEDs. 
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Fig. 8. — Images of edge-on disks at A = 0.814 fim for dust mixtures with p=3.5 (upper panels) and 
p=2.5 (lower panels), and 

O'max — 10 //m (left), CL max — 1mm (center), and CL max — 10cm (right). 
The model has = 0.046 M Q , and a central star with = 0.9 L Q and T* = 4000 K. Each panel 
is labeled with the corresponding p and a max . The images are calculated at A = 0.814 fim and 
convolved with the PSF of the HST PCI. The outermost contour level corresponds to 18.3 mag 
arcsec -2 , the peak contour level is 13.3 mag arcsec -2 , and the contour interval is 0.5 mag arcsec -2 
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Fig. 9. — The left panel is the image of HK Tau/c ( [Stapelfeldt et al. 1998j ) and the right panel is 
the image of a disk model with M = 3 X 1(T 8 M yr" 1 , a = 0.01, R d = 110 AU, M d = 0.065 M , 
M* = 0.5 Mq, i?* = 1.3 R , =3500 K, a max = 1 m, p = 3.5. The disk inclination angle i = 80.8° 
and the P. A. is 40°. The levels are the same than in Figure |8| and in the observed image. 
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Fig. 10. — Model disk flux at 1.3 mm vs. spectral index n between 0.769 mm and 1.3 mm (lines) 
compared with observed values (stars) for Taurus from BSCG and BS91. Lines connect models of 
the same mass: M d ~ 0.023 M (dashed line), M d ~ 0.046 M (solid line), and M d ~ 0.092 M Q 
(dotted line). Model values are calculated for maximum grain sizes a max = 10, 100, 120, 150, 300 
^m, 1 mm, 1cm and 10 cm, identified by dots along the lines. For reference, the dot corresponding 
to a max = 1 mm is encircled. Values for a max < 10 ^im are similar to those for a max = 10 (im. 
Upper panel: p = 3.5, lower panel: p = 2.5. Typical 1-er errors are shown in the upper left 
corner of the upper panel. The range of a covered by each sequence of models is a = 0.02-0.04 
(M d ~ 0.023 M ), a = 0.01-0.02 (M d ~ 0.046 M ), and a = 0.005-0.01 {M d ~ 0.092 M ). We also 
show the flux and spectral index for flat irradiated models with the same masses, and a max = 1 
mm (circled dot, at the lower left corner) and for well mixed models with ISM-dust (big triangles). 
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Fig. 11. — Model disk flux at 1.3 mm vs. spectral index n between 0.769 mm and 1.3 mm (lines) 
for different disk radius, inclination angle and mass accretion rate, (a) For the disk model with 
p =3.5, Md ~ 0.046 Mq and Rd = 100 AU, we change the disk radius (and consequently the 
mass) keeping M = 3 x 10 -8 M Q yr _1 and a = 0.01. The plotted values are: Rd =50 AU (dashed 
line), 100 AU (solid line) and 200 AU (dotted line). The values of imax and the observations are 
described in Figure |l(]. (b) For the disk model with p =3.5, Md ~ 0.046 M and Rd = 100 AU, we 
change the inclination angle, cosz = 0.25 (dashed line), 0.5 (solid line) and 0.9 (dotted line), (c) 
Non irradiated viscous disks. Each symbol corresponds to a different mass accretion rate (labeled 
by log M) and the lines connect models with the same mass, Md = 0.092 (dotted line), 0.046 (solid 
line) and 0.023 M (dashed line). For all these models a max = 1 mm, p = 3.5 and i = 0°. (d) The 
same models shown in panel (c), irradiated by the central star (assuming dust and gas well mixed). 
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